Asymmetric Explosion of Type la Supernovae as Seen from Near 
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ABSTRACT 

We present near-infrared spectra of late phase (> 200 d ) Type la supernovae (SNe la) taken at 
the Subaru telescope. The [Fe n] line of SN 2003hv shows a clear flat-topped feature, while that 
of SN 2005W show less prominent flatness. In addition, a large shift in their line center, varying 
from —3000 to 1000 (km s _1 ) with respect to the host galaxies, is seen. Such a shift suggests the 
occurrence of an off-center, non-spherical explosion in the central region, and provides important, 
new constraints on the explosion models of SNe la. 



Subject headings: infrared: stars 
2003hv, SN 2005W) 



supernovae: general — supernovae: individual (SN 2003du, SN 
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Introduction and Summary 



The brightness and approximate uniformity of 
Type la supernovae (SNe la) enable us to use 
them as reliable and distant standard candles, 
and provided evidence for the accelerating uni- 



verse ( Riess et al. 1998 : Perlmutter et al. 1999h . 
For precision cosmology, it is critical to under- 
stand the origin of diversity, as well as the ex- 
plosion mechanism of SNe la that is still open to 
debate (see, e.g., iHillebrandt fe Niemeverl l2000l 
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for a review) . Possible sources of diversity include 
the asymmetry of the explosion ([Wang et al. 2003) 



CWDs) dPiersanti et al. 2003: 


Saio & Nomoto! 


2004 


; Uenishi. Nomoto & Hachisu 


2003: Yoon & Langer 


2003 


) and /or turbulent behavior of the deflagra- 



tion flame (e.g., iGamezo et al. 2003; R opke et al 
120061) . Therefore, it is important to observation- 
ally investigate the distribution of the synthesized 
elements and the kinematical structure of the 
ejecta in order to constrain the explosion mod- 
els. 

In this respect, late-phase (~ lyr since the 
explosion) spectroscopy at near-infrared (NIR) 
wavelength provides important diagnostics. Be- 
cause the ejecta become optically thin in late 
phases, spectroscopy provides an unbiased, di- 
rect view of the innermost regions. Well isolated 
[Fe 11] emission lines at 1.257 and 1.644 /<m enable 
us to trace the distribution of the most important 
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isotopes synthesized in SNe la, i.e., Ni (which 
decays into 56 Co and then 56 Fe) and other iron iso- 
topes such as 54 Fe. NIR observations of the type la 
SN 1991T show a strongly peaked [Fe n] 1.644 jum 
line which appears symmetric wi th respect to the 



rest-frame of the host g alaxy (jSpyromilio et al 
119921 : iBowers et al.lll997ti . The NIR [Fe III line 



profile of SN 1998bu observed by lSpvromilio et all 
(|2004l ) appears to have a somewhat less centrally- 
peaked emission profile. 

This approach was highlighted by the discov- 
ery of a flat-top ped [Fe III 1-644 urn emission line 



in SN 2003du (|H5fhch et alJ 120041 ). Such a line 



profile requires a central hole in the kinematic dis- 
tribution of the radioactive ejecta. This sort of 
hollow radioactivity dis tribution is predicted in ID 



explosion models (e.g., Nomoto et al. 19941 for a 



review) in which the innermost regions are burned 
under high densities, p > 10 9 gcm~ 3 , and electron 
capture produces stable isotopes of 58 Ni, 54 Fe, and 
56 Fe, rather than radioactive 56 Ni. However, such 
a distribution is contrary to the predictions of 
state-of-the-art 3D deflagration simulations which 
predict large-scale turbulent mixing in the inner 
layers, and thus do not confine 56 Ni outside the 
central core. The IHoflich et all (|2004l ) NIR spec- 



trum of SN 2003du also indicate that the [Fe n] 
1.257, 1.644 /im emission lines show a blueshift of 
500—1, 000 km s _1 with respect to the host galaxy. 

However, since the observed sample size of late- 
phase NIR spectra is so small, it is not clear 
whether the flat-topped and blue-shifted lines 
are general properties of SNe la. Such obser- 
vations are extremely difficult because even rela- 
tively nearby and bright SNe la are faint in NIR 
at late phases (e.g., H = 20 ~ 21 for SN 2003du 
at ~ 1 year after the explosion). We therefore 
have conducted NIR observations of SNe la using 
the Subaru with OH airglow suppressor (Table 1) 
and obtained 2 more late-phase 1 — 2/im spectra 
of SN 2003hv and 2005W. 

Detailed explanation and discussion of the spec- 
tra will be presented elsewhere (Motohara et al., 
in prep.). We briefly report two interesting results 
in this letter. (1) The [Fe n] lines of SN 2003hv 
also show a clear flat-topped feature. (2) The [Fe 
n] lines of SN 2005W show less prominent flatness. 
(3) A large velocity shift of the line center, varying 
from —3000 to 1000 (km s _1 ) with respect to the 
host galaxies, exists. 



Considering the uniformity in optical bands 
around the peak luminosity, the existence of such 
inhomogeneity and asymmetry in SNe la is sur- 
prising, and provides important, new constraints 
on the explosion models. 

2. Observations 
2.1. SN 2003du 

SN 2003du was discovered by LOTOSS on 
2003 Apr. 22.4 in UGC 939 1 at about 15.9mag 
((Schwartz fe Holvorcem1l2003l). I t was confirmed 
to be a SN la bv iKotak et all (|2003l ) who re- 
ported that the optical spectrum resembles that 
of SN 2002bo about 2 weeks before maximum, and 
reached maximum light in the i?-band (is max = 
d ) on 2003 May 6.3 UT (JD 2452766.3 ± 0.5) 
(jAnupama. Sahu. fe ,Iosell2005h . 

As reported IHoflich et all (j2004 . JiJ-band 



spectroscopy of SN 2003du was carried out using 
OHS /CISCO (jlwamuro et al]200lHMotohara et al 
2002) at Subaru. For this work, we re-reduced the 
data taken on 2004 Feb 27.5 (+297 d ), which were 
processed in a standard procedure of flat fielding, 
sky subtraction, bad-pixel correction, and resid- 
ual sky subtraction. The flux was scaled using the 
ff-band photometry that was taken on 2004 April 
2.5, assuming that the change in the brightness 
is negligible at the late phase. Wavelength was 
calibrated using the standard pixel-wavelength re- 
lation of CISCO, of which the systematic error is 
estimated to be less than 0.5 pixels (< 3A). 

2.2. SN 2003hv 

SN 2003hv was discovered by LOTOSS on 2003 
Sep. 9.5 (UT) in NGC 1201 at about 12.5mag 
( Beutler fe Li 20031 ). and confirmed to be a SN la 
by the spectrum taken on 2003 Sep. 10.4 which 
resem bles that of SN 1994 D two days after maxi- 
mum (Dressier et al 1 l2003l) . We therefore assume 
maximum light in the £?-band to be JD2452891. 

Our Jii-band spectroscopy of SN 2003hv was 
carried out using OHS/CISCO on 2004 Oct 6 
(epoch +394 d ). The J-ff-band spectroscopy con- 
sists of 6 frames of 2000 sec exposure with an 0'.'5 
slit, providing a wavelength resolution of ~ 400. 
The A2 star SAO 169939 was observed after the 
target to correct the atmospheric and instrumen- 
tal absorption pattern. The data were reduced in 
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the same manner as that of SN 2003du. 



2.3. SN2005W 

SN 2005W was discovered on 2005 Feb. 1.4 



(UT) in NGC 691 at about 15.2mag (jNakano fe Li 
20051 ). and confirmed to be a SN la about a 
week before maximum on 2005 Feb 2.7 (UT) 
(jElias-Rosa et al.l 120051 ). The expansion velocity 
is measured to be ~ 11600 km s _1 from the Si n 
line. We assume maximum light in the B-band to 
be 2005 Feb 10 (JD 2453413). 

Our Jii-band spectroscopy of SN 2005W was 
carried out on 2005 Sep 12 (epoch +214 d ) by 
OHS/CISCO. It consists of 4 frames of 1000 sec 
exposure with an 0'.'5 slit. The A2 star HIP 20091 
was observed after the target to correct the atmo- 
spheric and instrumental absorption pattern. The 
data were reduced in the same manner as that of 
SN 2003du. 

3. Results 

The NIR spectra are shown in Figure [TJ It can 
be seen that all the observed SNe la exhibit strong 
[Fe n] 1.257 fxm and 1.644 /im lines. We discovered 
two important features. 

First, the line center of [Fe n], corrected for 
the redshift of the host galaxies, is not identi- 
cal for the three events. With respect to the 
rest wavelength of the line (1.644 /um), 2 SNe la 
(2003du, 2003hv) show a blueshift, while the other 
SN (2005W) shows a redshift. The corresponding 
velocity shifts relative to the hosts' rest frame are 
-2600 kms -1 , -2100 km s" 1 , and +1400 kms" 1 , 
for SNe 2003hv, 2003du, and 2005W, respectively. 

Such a velocity shift is confirmed by a mid- 
infrared(MIR) spectrum of SN 2003hv taken by 
the IRS/Spitzer Space Telescope at ~ +360 d af- 
ter maximum brightness (Fig. [2j Gerardy et al., 
in prep.). One of the strongest emission features 
is identified as the ground-state fine-structure line 
of [Co in] at 11.89 nm (a 4 F 9/2 - a 4 F 7/2 ). The 
well-isolated [Co in] line also shows a velocity shift 
which is consistent with that seen in the NIR [Fc 
n] 1.644 /im line, to within the noise level of the 
MIR feature. 

Secondly, we find that the observed [Fe n] lines 
show a large variety in their shape (Figure []}. 
The spectrum of SN 2003hv clearly shows a flat- 
topped boxy profile like that seen at much lower 



S/N in SN 2003du. Indeed, the higher S/N of 
the SN 2003hv spectrum places a much stronger 
constraint on the flatness of the core of the [Fe 
n] line. In contrast, the [Fe II] line profile of SN 
2005W shows no evidence of a central flat top. 

The boxy profile of [Fe n] lines seen in SN 
2003du and SN 2003hv suggests the absence of 
radioactive 56 Ni below ~ 3000 km s _1 (measured 
from the center of the 56 Ni distribution). We note 
that the observed epoch of SN 2003hv is more ad- 
vanced than SN 2003du. At a later epoch, the line 
profile should follow the 56 Ni distribution more 
closely (§4). Effects of possible line deformation 
due to electron and line scatterings will also ef- 
fectively vanish. The present result of SN 2003hv 
confirms the existence of the 56 Ni-empty hole more 
strongly than that of SN 2003du. 

4. Discussion 

In this section we discuss the impact of our find- 
ings for understanding the explosion mechanism. 
It is widely agreed that t he explosion of a SN la 
start s from a deflagration (jNomoto. Sugimoto. fe Neo 
19761 ). In the spherical deflagration models, elec- 
tron capture leads to the synthesis of 58 Ni, 54 Fe, 

and 56 Fe (not via 56 Ni dec ay) , thus creating an al- 

most 56 Ni-empty hole (e.g. jNomoto. Thielemann, fc Yokoi 
19841 ). In multi-D models, the ignition of the de- 
flagration may be off-center, producing a non- 
spher i cal burning region (IWunsch fc Wooslev 



2004 : 
cases. 



Plewa. Calder. fc Lamb! 120041) . In" 



some 



the deflagration to detonation transition 
(DDT) may oc cur (delayed detonation models: 
Khokhlovl[l99ll ) . The DDT may also take place 
non-spherically ( Livnelll999l ) even if the deflagra- 
tion does not have bulk kinematical offset. 

The flat-topped [Fe II] NIR lines in SNe 2003hv 
and 2003du indicate that the highest density re- 
gion occupied by neutron-rich Fe-peak isotopes 
is not mixed with the surro unding region wher e 
the dominant isotope is 56 Ni (jrloflich et al.ll2004th 
Figure [2] shows a spectrum for the spherical defla- 
gration explosion model W7 (Nomoto et al. 1984) 
at 400 d since the explosion, compared with that 
of SN 2003hv. This is calculated by solving trans- 
port of 7-rays produced in the decay chain 56 Ni 
-> 56 Co -> 56 Fe and iteratiyely solving NLTE 
rate equations (jMazzali et al. I l200lt iMaeda etajj 
20061 ). Positrons produced by the 56 Co-decay are 
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assumed to be trapped on the spot, since the mean 
free path o f positrons is exp ected to be small at 
this epoch (jMilne et alJl200lt l. 

The model shows that the flat-topped profile is 
consistent with the [Fe n] 1.644 fim emission for 
the reason as described above. The asymmetry in 
the profile, mildly peaking in the red, is due to the 
weak contributions of [Fe n] 1.664, 1.677 /im. 

The shift of [Fe n] lines suggests that the dis- 
tribution of 56 Ni (which decays to 56 Fe via 56 Co) 
produced at the explosion is asymmetric, showing 
the bulk kincmatical offset of > 2000 km s _1 with 
respect to the SN rest frame. This suggests that, 
unlike the spherical models, the carbon ignition 
may take place off-center. The flat-topped [Fe n] 
line profiles suggest that electron capture in the 
high density off-center deflagration region creates 
the neutron-rich hole as in ID models, at least in 
SNe 2003du and 2003hv. The neutron-rich hole is 
offset along with the bulk 56 Ni distribution. This 
indicates that the highest density burning in these 
SNe took place quite far away from the center of 
the progenitor star. 

The blueshifted and Hat-topped profile of [Fe n] 
can also be reproduced by symmetric and opaque 
dusty core, like [O i] AA6300, 6364 observe d in late 
phase Type II SN (|Elmhamdi et al.ll2003l ). If this 
is the case in the present results, MIR [Co in] line, 
which will be far less affected by dust extinction, 
may show much smaller blueshift than NIR [Fe 
ii]. However, the [Co in] line shows almost same 
blueshift as [Fe n] (Fig. [2]). Therefore, it is un- 
likely that the flat-topped profile is caused by the 
opaque dust core. 



Ad ding the result of SN 1991T (|Spyromilio et al 



119921: iBowers et al.lll997ft . 2 SNe (SN 2003du, SN 
2003hv) out of 4 SNe la show the boxy [Fe Ii] 
profile (Figured]), while SN1991T and SN 2005W 
do not. This could indicate that the distribution 
of 56 Ni in the innermost region of the SN la ex- 
plosion might differ from object to object. SNe 
la with a non-boxy line profile might experience 
more mixing than those with a boxy profile (e.g., 
SN 2003hv) in the innermost region. Another pos- 
sibility is that the density of the progenitor where 
the ignition takes place differs from object to ob- 
ject, and SN la with the peaked [Fe ii] explodes at 
a low density so as not to produce electron capture 
isotopes. 



Alternatively, the different line shape could be 
due to an age effect, i.e., due to a variation in the 
extent to which the energy deposition from the ra- 
dioactivity is kept local. At earlier epochs, 7-rays 
are the dominant heating source. Since the 7-ray 
penetration is not a local process, the innermost 
region, even without 56 Ni there, can be heated ef- 
fectively. At later epochs, contribution from the 
local positron energy input becomes larger, and 
the line profile should follow the 56 Ni distribution 
more closely. Therefore, we expect that the line 
shape could evolve from a peaked to a flat-topped 
profile. This may partially explain the fact that 
the most aged SN la 2003hv shows the clean flat- 
topped boxy [Fe 11] line. Indeed, this provides an 
observational test of our model. 

We suggest a future test to estimate the age 
effect and confirm our interpretation of the 56 Ni: 
a temporal series of NIR nebular spectra for an 
individual SN la. The shape of [Fe 11] line is ex- 
pected to start getting flat at ~ 250 d , and may 
become totally flat at ~ 500 d . Thereafter, it may 
start showing a peaked profile at ~ l,000 d again, 
depending on th e amount of positr ons that can es- 
cape the ejecta ( Milne et al.l 2001 ). Also, spectra 
at different wavelengths from optical to the MIR 
will be useful to investigate the distribution of dif- 
ferent ions and heating radioactive isotopes (e.g., 
Figure [2]) . Because the line profiles from differ- 
ent ions and from different energy levels are de- 
pendent on the ionization and thermal structure 
of the ejecta, further detailed theoretical study is 
necessary to make use of these observations effi- 
ciently. 
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Table 1 

Infrared spectroscopy observing log 



Name 


Host 


ZHoet 


Obs. Date (UT) 


Epoch" (d) 


AA(A) 


A/AA 


Exposure (s) 


Slit (") 


2003du 


UGC 9391 


0.006384 


2004 Feb 27.5 


+297.2 


11080-18040 


400 


4000 


0.5 


2003hv 


NGC 1201 


0.005604 


2004 Oct 06.5 


+394 


11080-18040 


400 


12000 


0.5 


2005W 


NGC 691 


0.008889 


2005 Sep 12.6 


+214 


11080-18040 


400 


4000 


0.5 



^Epochs arc with respect to the epoch of £?-band maximum light (tgrnax — 0). 
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\ Ml (Mm) 



Fig. 1. — NIR 1 — 2 /im spectra of late-phase SNe la, converted to the rest wavelength at the host galaxies. 
All the spectra are smoothed with 3 pixel boxcar filter. Vertical dashed lines show the position of [Fe n] 
1.257 fim and 1.644 /im lines. 
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Fig. 2. — [Fe Ii] 1.644 /jm (top) and [Co in] 11.89 fim (bottom) emission lines of SN 2003hv (solid). Over- 
lapped with the [Fe n] line is the ID explosion model W7 at 400 d since the explosion (dotted). The model 
flux is arbitrarily scaled to fit the observed flux, as well as shifted to the blue by 2600 km s _1 . 
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